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We study the influence of photons on the dynamics and the ground state of the atoms in a 
Bosonic Josephson junction inside an optical resonator. The system is engineered in such a way 
that the atomic tunneling can be tuned by changing the number of photons in the cavity. In this 
setup the cavity photons are a new means of control, which can be utilized both in inducing self¬ 
trapping solutions and in driving the crossover of the ground state from an atomic coherent state 
to a Schrodinger’s cat state. This is achieved, for suitable setup configurations, with interatomic 
interactions weaker than those required in the absence of cavity. This is corroborated by the study 
of the entanglement entropy. In the presence of a laser, this quantum indicator attains its maximum 
value (which marks the formation of the cat-like state and, at a semiclassical level, the onset of 
self-trapping) for attractions smaller than those of the bare junction. 


I. INTRODUCTION 


A Bose Einstein condensate confined in a double-well 
potential mimics in many ways the coherent dynamics of 
a superconducting Josephson junction HIE]. Therefore, 
it is often called a Bosonic Josephson junction (BJJ), 
whose coherent dynamics can be described by the nonlin¬ 
ear equations of a non-rigid pendulum [2|. Furthermore, 
the investigation of the BJJ dynamics allows us to study 
the formation of macroscopic coherent states mm and 
macroscopic Schrodinger’s cat states [iHinj. When the 
condensate is inside an optical resonator, the photon and 
atomic degrees of freedom are no longer independent m- 
there is a mutual back-action between the cavity and BJJ 
dynamics. The photon field acts as an optical potential 
on the atoms; at the same time, the atomic density serves 
as a refractive medium for the cavity photons. Note that 
the idea underlying this scheme is closely related to the 
models discussed, within the optical lattices context, in 
In an earlier paper it was shown, how the con¬ 
stant photon field of a laser - tightly focused to the center 
of the junction barrier - influences the tunneling proper¬ 
ties by effectively raising or lowering the barrier, or ul¬ 
timately allowing for a localized state inside the central 
well [in]. Later it was also investigated, how the dynami¬ 
cal nature of a cavity field - focused similarly to the center 
of the barrier - alters the BJJ dynamics in the semiclas¬ 
sical level [16]. The purpose of this paper is to study how 
the photons influence the exact quantum dynamics and 
the ground state of the junction, when the parameters 
of the system are chosen in such a way that the cavity 
field can be considered as fixed, i.e. in this bad cav¬ 
ity case the photon field is practically a laser field. The 
cavity photons alter the equations that describe the BJJ 
dynamics [16], and we investigate how they can induce 
self-trapping solutions. The study of the transition from 
the Josephson regime to the self-trapping regime in the 


semiclassical dynamics is complemented by the analysis 
of the ground state. The crossover from the atomic co¬ 
herent state to the Schrodinger’s cat state, can in fact be 
regarded as the quantum counterpart to the semiclassi¬ 
cal examination of the system. The ground state of the 
system is analyzed using three quantum indicators: the 
Fisher information, the coherence visibility and the en¬ 
tanglement entropy. This last estimator proves to play a 
central role in understanding when the transition takes 
place. 


II. THE MODEL 


Here we consider a hybrid system consisting of a BJJ 
and an optical cavity, whose field is interacting with the 
atoms. The scheme is identical to that of Ref. [T6] 
and is illustrated in Fig. [^ The BJJ in our system 
is made of an interacting Bose-Einstein condensate in 
a symmetric double-well potential. Its Hamiltonian in 
second-quantized form reads as 


Ha = J d^r 4'^(r) 


^\7^ + V{r) + ^¥{r)i>ir) 


4'(r), 


( 1 ) 

where ^'(r) annihilates an atom at position r = (x,y,z). 
Low energy atom-atom scattering is characterized hy g = 
ATTfi^a/m, where a is the s-wave scattering length, and m 
is the mass of an atom. The confining potential is 


V{r) = Vdw{x) z^) , (2) 


where Vdw{x) is the double-well potential in the x di¬ 
rection, while we assume tight harmonic confinement in 
the transverse direction, with trap frequency ujh- 

The optical cavity is a single-mode high-Q Fabry-Perot 
resonator, whose axis is orthogonal to the double-well di¬ 
rection as depicted in Fig. [^ We choose the y direction 
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FIG. 1. (Color online) The illustration of the setup. The 
bosonic Josephson Junction is created by magnetic or optical 
means along the x direction. A Fabry-Perot cavity is placed 
around the junction with an axis orthogonal to the junction. 
The resonator is operated on the TEMqo mode. 


for the cavity axis. The cavity has characteristic fre¬ 
quency ujCi it is pumped through one of its mirrors with 
a laser of frequency ujl and amplitude ry, and it is op¬ 
erated at the TEMqo mode. The mode function of the 
cavity is 

= cos(fc y) - - ^1/2^ ( 3 ) 

with k = wc/c is the wave number of the cavity mode, 
L is the distance between the mirrors and a is the width 
of the Gaussian profile in the {x, z) plane. 

The cavity dynamics, in the frame rotating with the 
frequency ujl of the laser drive, is governed by the Hamil¬ 
tonian 

He = — ihrj{a — a^). (4) 

Here Ac = ujl — wc is the cavity detuning, 77 is the 
pumping strength of the driving laser, d is the annihi¬ 
lation operator of a cavity photon and Nl = d^d is the 
photon number. 

The condensate atoms have an electronic ground state 
and an excited state, with transition frequency oja- We 
assume that the atomic detuning, = lol — is 
large enough to the excited state population to be neg¬ 
ligible all the time, and the atoms to behave as polariz¬ 
able scalar particles. Therefore, their interaction with the 
cavity field is described by dispersive photon scattering 
which shifts the cavity resonance frequency propor¬ 
tional to the atom number, and on the atomic degrees of 
freedom it creates an optical potential. This additional 
atom-light interaction term reads as 

Hal = KUod^d J d^r/^(r)4'^(r)^(r), (5) 


where the strength of the interaction is Uq = ilj^/A a, 
with the single-photon Rabi frequency. As the cavity 
field is focused symmetrically to the center of the double¬ 
well barrier, the optical potential creates a light shift 
symmetric in the wells, and lowers (raises) the height 
of the barrier for red (blue) detuned atoms. 

When the double-well barrier is sufficiently high, the 
lowest energy doublet is well separated from the rest of 
the spectrum. In this case the atomic dynamics can be 
constrained to the two-mode Fock space of the left and 
right valleys [2] and the atomic field operator can be ap¬ 
proximated as 

4 '(r) = (wi{x) bi + W2ix) b2j - ^1/21^ - ’ 


where bi and 62 are the annihilation operators of a bo¬ 
son in the Wannier-like states wi{x) and W2(x), cen¬ 
tered around the minima of the left and right valleys, 
respectively. The characteristic length of the strong 
harmonic confinement in the (y, z) plane is given by 
Ih = ^/h/{mitJH)■ In this limit, the Hamiltonian de¬ 
scribing the BJJ dynamics is evaluated by inserting the 
field operator ^ into the atomic Hamiltonian Q, 

i?Bjj = eN-J (blh + blbi)+^ {b\b\bih + blblhh) , 

( 7 ) 

The parameters e, J > 0 and U are the on-site energy 
of a single well, the tunneling amplitude and the on-site 
interaction energy, respectively. The values of these pa¬ 
rameters are written in terms of overlap integrals of the 
Wannier functions m- 

Similarly, the atom light interaction (§ evaluates to 


Hi = Nl 


W0N + W12 {b\b 2 + blbi^ 


( 8 ) 


The term with Wq is the symmetric AC-Stark shift, while 
the IT12 term is the consequence of the lowering of the 
barrier. This latter can be understood as a cavity assisted 
tunneling contribution. The value of the parameters Wq 
and IT12 are also expressed as overlap integrals [TB]. The 
full Hamiltonian of the system can be written as: 


H = iJBjj + Hc+Hi 

= (e + WoNl) N-{J- Wi2Nl) (blh + blh 
+ ^ (i\b\bibi + b\h\b2b‘^ - HAcNl - ihrj^d 


-f 




( 9 ) 


A more detailed derivation of the Hamiltonian can be 
found in Ref. [ 16 ] . In this configuration, Hq and IT12 can 
be either positive or negative since they are proportional 
to the atomic detuning A^ [TBj. If the atoms are red 
detuned, i.e. A^i < 0, both Wq , W12 < 0; while for 
blue detuning, i.e. A^ > 0 , one has Wq , IT12 > 0 . The 
geometry of the system has a remarkable advantage: the 
cavity photons can influence the tunneling of the atoms 
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between the two wells of the potential, as opposed to 
previous proposals of this type of system [nHTS! which 
considered only the effect of the cavity photons on the 
on-site energies of the BJJ. 


III. SEMICLASSICAL APPROXIMATION 

Let us assume that the system is in a full coherent 
state, in which both the atoms in the right and left 
well, and the cavity photons are described with coherent 
states: \FCS) = |/ 3 i) 0 |/ 32 ) ® la), where bj\l 3 j) = ( 3 j\/ 3 j) 
and a I a) = ala). It is convenient to write the eigenvalues 
of the atomic coherent state as: Pj = where 

Nj the average number of atoms in the j-th well and 9 j 
the corresponding phase, and similarly a = where 

Nl = is the average number of photons in the cav¬ 
ity and (j) the corresponding phase. The dynamics of the 
system can be described by a set of four ordinary differen¬ 
tial equations for the fractional imbalance of the atomic 
population of the wells z = ^ the atomic relative 

phase 9 — 92 — 9 i, and the photon variables ^ and 4 > 
[ 16 ]. The photon dynamics can be adiabatically elimi¬ 
nated when Sc = Aq — N{Wo + Wi2Vi — cos 6 *)/h is 
orders of magnitude larger than = (J — W\2^)lb? 
[ni. When the magnitude of the cavity detuning Aq is 
much larger than that of NW^jh and NWi2/h the sys¬ 
tem is described by: 

z = — 2 f'\/w^^sin 0 , ( 10 a) 

5 + cos 6 i) z, ( 10 b) 

where v = {J — W\2^)lh = (J — W\2'rf'jS' q) jh. When 
these conditions are satisfied, the photon held ampli¬ 
tude ^ can be considered to be constant in time and 
^ = r7/|Aq|. This set of equations is very similar to 
the Josephson equations for a bare Josephson junction. 
But here, because of the cavity photons, the tunneling 
amplitude J is substituted with the assisted tunneling 
amplitude J = J — Wi2^'^. In the case of the bare BJJ 
the solutions to these equations can be divided into two 
classes. One is characterized by oscillations around z = 0 
and therefore the time average of the population imbal¬ 
ance is zero (Josephson regime). The second class of 
solutions is characterized by an average imbalance in the 
population of the wells (Self-trapping regime). However, 
in our system, the cavity photons can be used to induce 
self trapping solutions for initial conditions and param¬ 
eters that would not allow them in the bare junction. 
When U > 0 and J > 0 , the self trapping occurs when 

m- 

|z2(0)- Vl-2"(O)cos0(O) > 1. (11) 

where A = UN/ 2 J and z( 0 ) and 0 ( 0 ) are the initial con¬ 
ditions to the equations. Since A depends on the number 




FIG. 2. (Color online) The highlighted area in the panels shows 
the values of z, U and 9 that allow self trapping. In the first column 
for a fixed 9 = 0 and in the second column for a fixed [/ = 12J/A. 
The difference between the two rows is the number of photons in 
the cavity. In the first one = 0 and J = J, in the second = 25 
and J = 0.25J. It is clear that, because of the cavity photons, 
the highlighted area gets wider and therefore the appearance of 
the self trapping is assisted. The other parameters are N = 1000, 
NW \2 = 30J. U is written in units of J. 


of cavity photons in our system we have one more 
relevant parameter that allows us to switch between dif¬ 
ferent regimes. By comparing the panels in Fig. it 
can be clearly seen what the effect of changing the num¬ 
ber of photons is when W12 > 0 . In particular we notice 
how increasing the number of photons widens the area 
that allows self trapping solutions. The photon-induced 
change of regime can be clearly seen from Fig. The 
left panel shows the time evolution of z in the Josephson 
regime, in the absence of photons. The right panel shows 
the time evolution when the cavity photons are present 
in the system: it can be clearly seen that the photons 
induce oscillations around a non-zero average value of z. 
The cavity photons can be therefore used as a means of 
inducing self trapped solutions that would not otherwise 
occur in the absence of the cavity. 


IV. EXACT DIAGONALIZATION 


We now want to show how the cavity photons influ¬ 
ence the ground state of the system. In order to do this, 
we describe the photons with a coherent state |a), but 
we describe the atoms in the wells in a purely quantum 
way. To this end we use Fock states as basis vectors 
in the atomic Hilbert space. Furthermore, we assume 
that the photon field relaxes quickly and this allows us 
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FIG. 3. (Color online) Time evolution of the variable z. Time is 
measured in units of h/ J. The first panel represents the bare junc¬ 
tion in the absence of the cavity. The second one shows the junction 
inside of the cavity. The solid lines represent the numerical solution 
of the system where the photon variables are constant. The 

dashed line represents the solution of the system where the photon 
dynamics has not been adiabatically eliminated. When there are 
no photons, the system is in the Josephson regime, while when the 
photon number is increased the system is in the self trapped regime. 
The parameters for the second panel are: hAc = 300J, WqN = 4J, 
W 12 N = 3J, UN = 12J, N = 1000, hr) = 3000J. The initial con¬ 
ditions for the first panel are for ( 2 , 9) = (0.7, 0). The initial condi¬ 
tions for the second panel are {z,9,^,<j>) = (0.7,0,10.17,7r/2). The 
inset plot represents the frequency uj of the small oscillations of 2 
written as a function of the number of cavity photons . ui is writ¬ 
ten in units of JI h. The range of is given by the condition for the 
system to show small oscillations. The other relevant parameters 
for the inset are N = 1000, UN = 0.012J and W 12 N = 0.03J. 


to consider the photon coherent state to be a parameter 
on which the atomic ground state depends. The ground 
state of the atomic held for a given photon held coher¬ 
ent state |a) can be written as: \GS)a = 
where \n) = \N — n,n) is an element of the Fock basis 
{|n), n = 0,1,..., N} and N is the total number of atoms. 
The element of the Fock basis \n) describes the state in 
which n atoms are in the right well and iV — n in the left 
well. The coefficients c„(a) clearly depend on the photon 
held. So the ground state of the system can be written 
as IGS”) = |a) ® |GS')q, which is the tensor product of 
the coherent state of the photon held and of the atomic 
ground state for that held strength. In order to calculate 
the ground state coefficients c„(a) the matrix elements 
of the Hamiltonian that needs to be diagonalized are: 


= - (J - IFmlaH 

-I- y/N - nv^5m,n+i) + ^[(-^ - n)(A^ - n - 1) 

-I- n(n - l)]Jm,n, 

( 12 ) 


If we assume that the photon held relaxes very fast, then 
ih = 0 leads to a relation between a and the ground 
state coefficients c„(a) that needs to be satished in order 
for the computation of the ground state to be consis¬ 
tent. This makes the amplitude of the photon held |a) 
dependent on the parameters that appear in the Hamil¬ 
tonian and on the coefficients c„(a) as well. The anal¬ 
ysis of the ground state can however be simplihed when 
I Ac I ^ \Wo\N/h and |Ac| S> \Wi 2 \N/h. When such 
conditions are satished a depends weakly on the coeffi¬ 
cients Cn{ot). Therefore, for any starting number of pho- 
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FIG. 4. (Color online) Coefficients |cn(o:)P of the ground state 
as a function of n/N for N = 1000. In panel (a) W \2 = —0.03J 
and in panel (b) W \2 = 0.03J. The on-site interaction is nega¬ 
tive f/ < 0. The solid line represents the system in the absence 
of cavity photons, \a\^ = 0, while the dashed line represents the 
system when |op = 1. The effect of the cavity photons in (a) is 
to increase the coherence of the ground state for a given negative 
on-site interaction I/, and thus delay the transition to the regime 
in which a valley appears in the middle of the distribution of the 
coefficients |cn(o:)p. The effect of the cavity photons in (b) is the 
opposite. The coherence of the ground state decreases for a given 
negative on-site interaction [/, and thus the photons assist the tran¬ 
sition to the regime in which a valley appears in the middle of the 
distribution of the coefficients |cn(a)|^. In both cases the ground 
state is showed for four values of U. The coefficients |cti(o)P are 
adimensional and normalized so that \cn{(y.)\‘^ = 1. n/N is 

adimensional. 


tons lap the above procedure gives as result the starting 
|ap. The number of photons in the system can thus 
be modified by changing the values of rj and Ac. Under 
this approximation we only need to study how the ground 
states of the Hamiltonian with the assisted tunneling is 
influenced by the presence of the photons. In the ab¬ 
sence of the cavity, and thus |ap = 0 the ground state of 
the BJJ can show three limit behaviors depending on the 
value of the ratio U/J [51]. When U/J = 0 the ground 
state is an atomic coherent state (ACS), where the dis¬ 
tribution of the |c„(a)p has a gaussian profile. For large 
values \U/J\ the ground state can be found in two config¬ 
urations, depending on the sign of U/J. When U/J>0 
the ground state is a separable Fock state, where half of 
the atoms are in the right well and half in the left one. 
When U/J<0 the ground state is a ’’cat state”, which 
is the linear combination of the state with all particles 
in the left or in the right well. Changing the interaction 
strength U/J leads to a continuous transitions between 
these three kinds of state. 

The correlation properties of the ground state can be 
analyzed by using the Fisher information F (53] and 
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FIG. 5. (Color online) Fisher information F, coherence visibility 
OLy and entanglement entropy S, plotted as a function of the on¬ 
site interaction 17, which is written in units of J, for N = 1000. 
In panel (a) VF 12 = —0.03J and in panel (b) VF 12 = 0.03J. The 
solid line represents the system when |ap = 0, while the dashed 
line represents the system when \a\^ = 1. The on-site interaction is 
attractive U < 0. The plots of the quantum indicators complement 
those of Fig. In the presence of photons F is smaller when 
W 12 < 0 and larger when W 12 > 0, compared the the solid line, 
since F is related to the distribution width of the Cn(o:). ay is larger 
for W 12 < 0 and smaller for VF 12 > 0 compared to the solid line 
as the photons increase and decrease the coherence respectively. 
As for S, whose maximum represents the value of U at which the 
distribution of the |cn(»)|^ starts to show a valley m, comes for 
a larger \U\ for W 12 < 0 and for a smaller \U\ for VF 12 < 0. F,ay 
and S are dimensionless. 


the coherence visibility q;„ [21] • The definition for these 
two estimators can be found in m- In particular, F gives 
the width of the distribution of |c„(a)p and charac¬ 
terizes the degree of coherence between the two wells. We 
are interested in characterizing also the genuine quantum 
correlations between the atoms in two wells pertaining 
to the ground state. Here we focus on the atom part 
only and address this issue from the bi-partition per¬ 
spective with the two wells playing the role of the two 
partitions. Following the same path as in [H US] we 
thus calculate the entanglement entropy S [21] that re¬ 
sults in S' = -Yln=o |cn(a)plog 2 |c„(a)p. Remarkably, 
S, plotted as a function of U, shows a maximum at the 
onset of the transition to the ’’cat-like” state regime in 
which the distribution of the |c„(a)p begins to exhibit 
a valley. The Hamiltonian we need to diagonalize in Eq. 
(121 differs from the Hamiltonian of a Bosonic Josephson 
Junction only for the value of the tunneling amplitude, 
which we will call J = J — kFi 2 |ap. The magnitude and 
sign of J depend on both W 12 and on the photon number 
lap. 


When W 12 < 0 the assisted tunneling amplitude J is 
always positive and its magnitude gets larger as the num¬ 
ber of photons in the system |ap increases. Let us choose 


W 12 = —0.03J and plot together the ground states for 
various values of ?7 for |ap = 1 and therefore J = 1.03J. 
We always plot the ground state of the bare junction, 
with jap = 0, to use as a reference. In panel (a) of 
Fig. [4] and of Fig. [^we can see what happens when the 
on-site interaction U is negative. When the interaction 
is positive, the effect of the photons is still to delay the 
transition to the separable Fock state. 


When W 12 > 0 the assisted tunneling amplitude J 
gets smaller when increasing the number of photons in 
the cavity. This can lead to two main consequences. On 
the one hand, if the number of photons is sufficiently 
small, J maintains its positive sign but its magnitude 
gets smaller and 0 < J/J < I. On the other hand if the 
number of photons is large enough the assisted tunnel¬ 
ing can become negative J/J < 0. When 0 < J/J < 1 
we can plot the ground states of the system, for differ¬ 
ent values of negative U, and for |ap = 1, and therefore 
J = 0.97J. We can see from panel (b) of Fig. and 
of Fig. that the effect of the photons in this case is 
opposite to the effect when IF 12 is negative. When the 
interaction is positive, the effect of the photons is still 
to assist the transition to the separable Fock state. It is 
worth mentioning that this result is actually more inter¬ 
esting then the sole rescaling of the ratio U/J. Focusing 
our attention on the attractive on-site interaction and 
on the emergence of the ’’cat-like” state regime we can 
notice a fundamental difference. For a fixed number of 
atoms in the system, in the bare Josephson junction, the 
transition to the ’’cat-like” state regime happens for a 
definite value of U/J. However, in the complete system, 
with a positive W 12 , by fine tuning the number of cav¬ 
ity photons the transition to the ’’cat-like” state regime 
can happen for a given magnitude of negative U, without 
changing the number of atoms. This means that the pho¬ 
tons in the system act as a new knob through which the 
transition between different regimes can be manipulated. 
The change of the sign of J is significant when analyzing 
the Hamiltonian of the system. In the case of the bare 
Josephson junction J can always be taken to be positive. 
However, by numerically calculating the ground state of 
the system we found out that the sign of J bears no rel¬ 
evance to the coefficients |c„(a)p and to the value of the 
quantum indicators. Therefore, choosing values of |q;^| 
such that the values of the assisted tunneling amplitude 
J is the same as the ones we have just studied, but with 
opposite signs, leads to the same results. It is interest¬ 
ing to notice this analogy between the case with positive 
J and negative J. One would at first think that, since 
the relevant parameter in the analysis is the ratio U/J, 
when J is negative, the same results could be obtained 
by changing the sign of U as well, however this is not the 
case. 
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V. SUMMARY 

To summarize, we have shown that the cavity photons 
can be used to induce self-trapping solutions in the dy¬ 
namics of the Josephson junction. Moreover we have de¬ 
scribed how the transition from the atomic coherent state 
to a ’’cat state” can be delayed or expedited, depend¬ 
ing on the photon-assisted tunneling amplitude. This 
analysis lays the first stone of the quantum study of this 
system. Possible extensions to this work could be the 
study of the ground state when the |Ac| 3> \Wi)\N/h 
and I Ac I |lTi 2 |A^/^ conditions no longer hold. This 
would lead to the same ground states we have already 
found, but the photon number of the cavity would de¬ 
pend on the other parameters of the system. Moreover 
one could try to better understand the transition between 
the atomic coherent state and the cat state, by comput¬ 
ing the relevant critical indices for F and 

The setup considered here can be implemented directly 
due to the current feasibility of trapping ultracold bosons 
in double-well potentials and to fabricate high- 

finesse optical cavities supporting quantum degenerate 
bosonic clouds [30H32j . The experiment deriving from 
the suitable combination of these two parts would pro¬ 
vide a powerful tool for observing the Josephson-self- 


trapping crossover, an issue addressed with bare BJJs 
m and creating entangled states of matter. In this con¬ 
text, experimentalists could have the concrete possibility 
to receive feedback on the shift of the crossover above by 
using the number of cavity photons as a knob. The setup 
considered in Ref. [31] is extremely appealing, where 
the Fabry-Perot cavity is assembled on top of an atom 
chip, where wires suitable for the double-well potential 
can be easily placed too. On such a device one has to 
fine tune experimental parameters that mainly consist 
in the proper matching of the double-well dimensions (a 
barrier width of a few microns) with those of the waist 
of the cavity (also in the range of a few microns). 
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